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Abstract 

We consider planar quadrangulations with three marked vertices and discuss 
the geometry of triangles made of three geodesic paths joining them. We also 
study the geometry of minimal separating loops, i.e. paths of minimal length 
among all closed paths passing by one of the three vertices and separating 
the two others in the quadrangulation. We concentrate on the universal scal- 
ing limit of large quadrangulations, also known as the Brownian map, where 
pairs of geodesic paths or minimal separating loops have common parts of 
non-zero macroscopic length. This is the phenomenon of confluence, which 
distinguishes the geometry of random quadrangulations from that of smooth 
surfaces. We characterize the universal probability distribution for the lengths 
of these common parts. 



1. Introduction 

Understanding the geometry of large random quadrangulations is a fundamental 
issue relating combinatorics, probability theory and statistical physics. Indeed random 
quadrangulations, or more generally random maps, provide natural discrete models for 
random surfaces, for instance in the context of two-dimensional quantum gravity 
and may mathematically be viewed as metric spaces endowed with the graph distance. 
In the same way that discrete random walks converge to the Brownian motion in a 
suitable scaling limit, it is expected that random planar quadrangulations converge to 
the so-called Brownian map in the scaling limit where the size of the quadrangu- 
lation becomes large jointly with the fourth power of the scale at which distances are 
measured. This Brownian map is moreover expected to be the universal scaling limit 
for many models of planar maps such as random planar triangulations or more gen- 
erally maps with arbitrary bounded face degrees or even maps coupled to non-critical 
statistical models. It can be constructed as a random metric space and has been shown 
to be homeomorphic to the two-dimensional sphere [^0]. 



A number of local properties of the Brownian map can be derived from a detailed 
analysis of discrete maps. In this spirit, the simplest observable in the Brownian map 
is the distance between two points. The statistics of this distance is characterized by 
the so-called two-point function and was obtained in Ref. via scaling arguments for 
large triangulations, and in Ref. P] via an exact computation of the discrete two-point 
function for planar quadrangulations. A related quantity is the radius, whose law was 
studied in Refs. [|TO| , p!T[] . The question of estimating the number of geodesies (i.e. paths 



of shortest length) between two points was addressed later ||T^ and it was found that 
for typical points, all geodesies coalesce into a unique macroscopic geodesic path in the 



scaling limit [13,14 
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Fig. 1: A schematic picture of the phenomenon of confluence for the 
geometry of triangles (a) and separating loops (b) in the scaling limit of 
large maps. In (a), the three geodesies (represented as thick blue lines) 
linking the three points Vi, V2 and ^3 have common parts of macroscopic 
length. The triangle is therefore characterized by six lengths (as indicated 
by double arrows) and by the area of the two domains delimited by its 
open part. In (b), a minimal separating loop passing by ^3 and separating 
Vi from V2 also has a common part of macroscopic length and is therefore 
characterized by two lengths (as indicated by double arrows) and by the 
area of the two domains delimited by its open part. 



Properties involving three points on the map give a much richer geometric informa- 
tion. For instance, we may consider the "triangle" made by the three geodesies between 
these points. In a previous paper [|15|, the authors have computed the joint probability 



distribution for the pairwise distances between three uniformly chosen random vertices 
in a random quadrangulation. In the scaling limit, this yields the so-called three-point 
function of the Brownian map, which can interpreted as the joint law for the three side 
lengths of the triangle. On the other hand, it was recognized by Le Gall that geodesies 
exhibit a phenomenon of confluence ||14|| . In our setting, this means that any two sides of 



the triangle merge before reaching their common endpoint, hence have a common part 
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of non-zero macroscopic length. This is quite unhke smooth surfaces where two sides of 
a triangle only meet at their endpoint. Thus a full characterization of the geometry of 
triangles involves six lengths, which are those of the three segments proper to each side 
and of the three segments common to two sides (see Fig. l-(a)), as well two areas for 
the two domains in the map delimited by the triangle. 

Beyond triangles, another interesting geometric construction involving three points 
is what we call a minimal separating loop, defined as follows: given three distinguished 
points, say fi, V2 and ^3, we define a separating loop as a closed path passing through 
f3 and separating vi from V2, in the sense that any path from vi to V2 necessarily 
intersects it. A minimal separating loop is such a separating loop with minimal length. 
We expect the minimal separating loop to be unique at a macroscopic level, and to 
have a finite macroscopic length (note that, if we relax the condition that the loop 
passes through ^3 or that it separates vi from ^2, then clearly we can find loops of 
arbitrary small length). Moreover, its two halves are geodesic paths and we again 
expect a phenomenon of confluence, namely the two halves share a macroscopic common 
segment (see Fig. l-(b)). The characterization of the geometry of minimal separating 
loops involves therefore the lengths of its common and "open" parts, as well as the areas 
of the two domains delimited by the loop. 

In this paper, we derive the probability distributions for the above parameters char- 
acterizing triangles and loops when the three points are chosen uniformly at random. 
This is done by explicit computations of the discrete counterparts of these distribu- 
tions in the framework of planar quadrangulations, using the methodology developed 
in Ref. and based on the Schaeffer and Miermont bijections between quad- 
rangulations and well-labeled maps. 

The paper is organized as follows: in Section 2, we give a precise definition of 
minimal separating loops in triply-pointed planar quadrangulations and compute the 
generating function for such quadrangulations with a prescribed value for the loop length. 
To this end, we provide in Section 2.1 two alternative bijections based on the Schaeffer 
and Miermont constructions relating the desired class of triply-pointed quadrangula- 
tions with suitable classes of well-labeled trees or maps. In Section 2.2, we calculate 
their generating functions by expressing them in terms of basic building blocks already 
computed in Ref. [|T3[. Section 2.3 is devoted to the analysis of the scaling limit, with a 



particular emphasis on the universal probability law for the length of the minimal sepa- 
rating loop, as well as its correlation with the distances between the marked vertices. In 
Section 3, we turn to the phenomenon of confluence, which we analyze by a refinement 
of the above enumeration. In Section 3.1, we give the probability law for the length 
of the part common to two geodesies leading to the same vertex. We then investigate 
the phenomenon of confluence for minimal separating loops in Section 3.2 where we 
derive the probability distribution for the parameters characterizing the geometry of 
these loops. Section 4 is devoted to the geometry of triangles. There we revisit the 
bijection of Ref. and solve a refined enumeration problem in order to keep track 
of the six lengths characterizing the triangle. We deduce their joint law in the scaling 
limit, and provide explicit expressions for a number of marginal laws. We discuss our 
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(a) 



(b) 



Fig. 2: A quadrangulation with tliree marlced vertices vi, V2 and ^3, 
represented as a traffic network, i.e. a ribbon grapfi witfi roundabouts. In 

(a) , we sliow (red tliick lines) a particular separating loop of length 6. In 

(b) , the indicated separating loop is minimal, i.e. has a minimal length 
(here 4) among the loops passing through ^3 and separating vi from V2- 

results and conclude in Section 5. 

2. Minimal separating loops 

Consider a quadrangulation of the sphere, i.e. a planar map whose faces all have 
degree four, equipped with three marked distinct vertices f 1, V2 and ^3. As customary 
for orientable maps, we may represent the map as a ribbon graph by splitting each edge 
of the quadrangulation into two oriented half-edges (with opposite orientations) so that 
half-edges are oriented clockwise around each face (see Fig. 2). It is also convenient 
to place a small counterclockwise oriented roundabout around each vertex so that the 
map looks like a traffic network. We can then consider (oriented) paths on this traffic 
network, and in particular loops made of a closed non- intersecting circuit starting from 
and returning back to the marked vertex ^3. Any such loop separates the sphere into 
two simply connected domains. Note that any vertex along the loop naturally belongs 
to exactly one of these domains by following the roundabout convention. The circuit is 
called a separating loop if the marked vertices vi and V2 do not lie in the same domain 
(see Fig. 2 for an illustration). The length of a circuit is the number of half-edges it 
passes through. A minimal separating loop is a separating loop of minimal length. 

Clearly the length /123 of a minimal separating loop is strictly positive and, from the 
bipartite nature of planar quadrangulations, it is even. Also, if we call di3 (respectively 
d23) the graph distance from vi (respectively V2) to ^3, following a geodesic path back 
and forth from ^3 to the closest vertex vi or V2 forms a separating loop of length 
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2 min((ii3, ^23)5 therefore: 

Z123 < 2min(rfi3,(i23) 



(2.1) 



The purpose of the next sections is to enumerate triply-pointed quadrangulations whose 
three marked vertices have prescribed values of (iis, ^23 and /123. 

An alternative definition of separating loops, mentioned in the introduction, con- 
sists in taking arbitrary (possibly self-intersecting) closed paths passing through ^3 and 
such that any path from vi to f 2 necessarily intersects them. This gives rise to a broader 
set of minimal separating loops but does not affect the minimal length since any such 
minimal separating loop can be transformed into a non self-intersecting circuit of the 
same length by "undoing" the crossings. 

2.1. Combinatorics 

Approach via the Schaejfer bijection 




Fig. 3: The quadrangulation of Fig. 2 with a marked origin (corresponding 
to vs in Fig. 2) and its coding (a) by a well-labeled tree (blue thick lines). 
The quadrangulation is recovered from the well-labeled tree by connecting 
each corner to its successor (dashed red arrows in (b)). 



It is well-known [O] that any planar quadrangulation with n faces and a marked 



origin vertex is in one-to-one correspondence with a well-labeled tree with n edges and 
with minimal label 1. Here we define a well-labeled tree as a plane tree with vertices 
carrying integer labels £ satisfying 

\i{v) — i{v') \ < 1 if f and v' are adjacent in the tree . (2.2) 



As shown by Schaeffer [16], this tree can be drawn directly on the quadrangulation 



by applying local rules which associate to each face of the quadrangulation an edge 
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of the tree (see Fig. 3). The tree spans aU vertices of the quadr angulation except the 
origin, and the label of each vertex is nothing but its graph distance to the origin in 
the quadrangulation. Conversely, to recover the quadrangulation from the well-labeled 
tree, we draw non-crossing arches connecting every corner of the tree to its successor. 
Recall that a corner is the sector between two consecutive edges around a vertex, and 
the successor of a corner with label £ > 1 is the first corner with label t — 1 encountered 
after it clockwise along the contour of the tree, while all corners with label 1 have the 
same successor which is an extra vertex added in the external face (see Fig. 3-(b)). 
The arches form the edges of the quadrangulation and the added vertex is the origin. 
Note that the chain of successors of a given corner (i.e. its successor, the successor of 
its successor, and so on until the origin) provides a geodesic path from the associated 
vertex to the origin. 




Fig. 4: A well-labeled tree with two marked vertices vi and V2- The edges 
of the branch from vi to f 2 are represented as magenta thick lines and the 
other edges as light-blue thin solid lines. The vertex is the origin added 
in the external face. We consider a vertex of minimal label (here 2) on the 
branch from vi to V2 and represent the chains of successors (dashed red 
arrows) starting from two of its corners, one on each side of the branch. 
These form a minimal loop separating vi from V2 and passing through ^3. 

In the case of a triply-pointed quadrangulation, we can take f 3 as the origin vertex 
and we end up with a well-labeled tree with two marked vertices vi and V2 carrying 
labels £(fi) = (ii3 and £(^2) = (i23- Let us now explain how the quantity /123 can be 
read off the tree. Within the tree, there is a unique branch connecting vi to V2 (see 
Fig. 4). Any loop separating vi from V2 in the quadrangulation must intersect this 
branch at some vertex v. Decomposing the loop into a first part from to v and a 
second part from v back to V2,, both parts have length larger that the distance l{v) from 
V to and we find that the length of the loop is larger than 2l{v), hence larger than 
2u, where u is the minimal label encountered along the branch from vi to V2- This holds 
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in particular for minimal separating loops, and we therefore have /123 > 2^. Conversely, 
a separating loop of length 2u is obtained by considering a vertex with minimal label 
u on the branch, picking two corners on opposite sides of the branch and considering 
the chain of successors of these two corners which are both paths to f 3 of length u (see 
Fig. 4). This implies /123 < 2u and therefore /123 = 2u. More generally, any minimal 
separating loop crosses the branch at a vertex with minimal label u, hence it is made 
of two geodesic paths of the same length u joining the origin to that vertex which they 
reach from both sides of the branch. 




i{v2) = t 



i{vi) = S 



mm£[v) = 1 



Fig. 5: The well-labeled tree coding a triply-pointed quadrangulation with 
prescribed values of dis, and /123. It has two marked vertices vi and V2 
with respective labels s = dis — /123/2 and t = ^23 — ^123/2. The minimal 
label on the branch between vi and f 2 is and the global minimal label is 
1-M = 1-/123/2. 



For consistency with the alternative approach described below, we decide to shift 
all labels on the well-labeled tree by —u so that the minimal label on the branch from 
f 1 to V2 becomes 0. The minimal label in the whole tree is now 1 — -u, while vi and 
V2 receive respective non-negative labels s = di3 — u and t = ^23 — u (see Fig. 5 for 
an illustration). To conclude, triply-pointed quadrangulations with prescribed values of 
di3, d23 and ^123 are in one-to-one correspondence with well-labeled trees having two 
marked vertices labeled s = di3 — /123/2 and t = ^23 — ^123/2, such that the minimal 
label on the branch joining these two vertices is and the global minimal label in the 
tree is 1 — w = 1 — /123/2. 



Approach via the Miermont bijection 

An alternative approach is based on a bijection by Miermont [13| generalizing the 
Schaeffer bijection to multiply-pointed planar quadrangulations. More precisely, the 
Miermont bijection acts on a quadrangulation equipped with, say p marked vertices 
f 1, f2, . . . fp, called sources and p integers ri, r2, . . . Tp called delays, satisfying the con- 
ditions: 

\Ti-Tj\<dij, l<i^j<p, ^ ^ 

(2.3) 

Ti — Tj + dij IS even, I < i, j < p , 
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(a) 



(b) 



Fig. 6: The quadrangulation of Fig. 2 witli ttiree marlced vertices vi, V2, 
V3, and its coding (a) by a well-Iabeled map (blue thick lines) using the 
Miermont bijection with particular delays ti = T2 = —1 and rs = —2. The 
quadrangulation is recovered from the well-labeled map by connecting each 
corner to its successor (dashed red arrows in (b)). 

where dij is the graph distance between Vi and Vj . It results into a planar map with p 
faces that is well-labeled, i.e. its vertices carry integer labels £ satisfying 

\i{v) — i{v')\ < 1 if f and v' are adjacent in the map. (2-4) 

Again, this map can be drawn directly on the quadrangulation by applying local rules 
which associate to each face of the quadrangulation an edge of the map (see Fig. 6). 
The map spans all vertices of the quadrangulation except the p sources and the label of 
a vertex v is given by 

£{v) = min d{v,Vj)+Tj (2-5) 

where d{v, Vj) is the graph distance from v to the source Vj in the quadrangulation. Each 
face of the well-labeled map encloses exactly one source of the quadrangulation and we 
call the faces /i, /2, • ■ ■ , /p accordingly. We furthermore have the property that, for any 
vertex v incident to fi, the minimum in ( pT5| ) is attained for j = z, i.e. d{v, Vi) = i{v) —Ti. 
In particular, the minimal label among vertices incident to fi is + 1, corresponding 
to nearest neighbors of Vi. 

Conversely, to recover the quadrangulation from the well-labeled map, we add inside 
each face fi an extra vertex with label Ti where 

r, = min e{v) - 1 , (2.6) 

V incident to fi 

and each corner with label i inside fi is connected by an arch to its successor, which is the 
first corner with label i — 1 encountered counterclockwise inside the face (corresponding 
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for the external face to the clockwise orientation around the map). The arches form the 
edges of the quadr angulation and the added vertices are the sources (see Fig. 6-(b)). 

Let us now see how to use the Miermont bijection to address the specific question of 
three marked vertices with prescribed values of dis, and Z123. As in Ref. [jl5|, the idea 



is to supplement the Miermont bijection (here with p = 3 sources) by a particular choice 
of delays related to di3, d23 and Z123. This particular choice will restrict the topology 
of the resulting well-labeled maps with 3 faces, and induce extra conditions on labels. 
More precisely, from the inequality (|2.1D , we may use the following parametrization: 



dis = s + u , 

d23 = t + u, (2.7) 

ll23 = 2u , 

with s, t, u non-negative integers, and moreover w 7^ 0. Our particular choice of delays 
is: 

Ti = -S = /123/2 - di3 , 

T2 = -t = /123/2 - ^23 , (2.8) 
T3 = -U = -/123/2 . 



Note that this particular choice fulfills the general condition ( |2.3| ) except when we have 
the equality /123 = 2min((ii3, (^23), i-e. when s or t vanishes. This particular case must 
be treated separately, as will be explained below. 

Assuming s and t strictly positive, a close look at the properties resulting from the 
choice of delays ( |2.8| ) in the Miermont bijection shows that the resulting well-labeled map 
is necessarily of the type displayed in Fig. 7. In particular, we find that any minimal 
separating loop must remain inside the face /a, hence the faces fi and f2 cannot be 
adjacent, i.e. cannot be incident to a common edge. The map can be viewed as made 
of a skeleton map (thick lines and big dots in Fig. 7) to which trees are attached. The 
skeleton is necessarily made of two cycles ci and C2 , which form respectively the frontier 
between /i and /s, and between /2 and /s, together with a bridge b connecting ci to C2, 
and whose edges are only incident to /a. Moreover, the labels must satisfy the following 
constraints (see Fig. 7): 

min i{v) = 1 — s , min i{v) = 1 — t , min i{v) = 1 — u , 

V incident v incident v incident 

to /i to /2 to /3 ^2 

min i{v) = , min i{v) = , min i{v) = . 

V on ci V on C2 v on b 

The first three constraints are general consequences of the Miermont bijection and 
rephrase the general condition ( p.6|) , while the last three constraints result from our 
particular choice of delays, and can be obtained by arguments similar to those pre- 



sented in Ref. [15|. More precisely, the constraint on ci (respectively C2) ensures that 



the distance between vi and V3 (respectively V2 and V3) is s + u (respectively t + u), 
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Fig. 7: The well-labeled map with three faces coding a triply-pointed 
quadr angulation with prescribed values of dis, and /123, in the case 
/i23 < min((ii3, (i23)- The faces fi and /2 are not adjacent, and their 
frontiers with the face fs form two cycles ci and C2, connected by a bridge 
b, whose edges are adjacent to fs only. The minimal label for vertices 
incident to /i (respectively /2 and fs) is 1 — s = 1— ^13-1-/123/2 (respectively 
1 — t = 1 — d23 + /123/2 and 1 — u = 1 — /123/2). The minimal label on the 
cycle ci is 0, as well as that on the cycle C2 and that on the bridge b. 



while the constraint on b ensures that the length of a minimal separating loop is 2u. 
Note that the bridge b can be reduced to a single vertex, necessarily with label 0. 

When s — and t > 0, we apply the Miermont bijection with p = 2 sources only, 
namely V2 and ^3, and delays T2 = —t, = —u. We obtain a well-labeled map with 
two faces of the type illustrated in Fig. 8- (a). In particular, vi is necessarily incident 
to /3 and has label 0, and is connected to the frontier between /2 and by a bridge 
having non-negative labels only. This can be seen as a degenerate version of the generic 
case displayed in Fig. 7, where the face /i is shrunk into a single vertex. We have a 
symmetric picture when s > and t = 0. Finally, if s = t = 0, we apply the Miermont 
bijection with p = 1 source only (equivalent to the Schaeffer bijection), namely f 3, and 
delay t^, = —u. We then obtain a well-labeled tree on which the vertices vi and V2 
have label and the branch connecting them has non-negative labels. Again this is a 
degenerate case of the generic situation in which both /i and /2 degenerate to single 
vertices. 

To conclude, triply-pointed quadrangulations with prescribed values of di2, (ii3 
and /i23 are in one-to-one correspondence with well-labeled maps of the generic type 
displayed in Fig. 7, or of its degenerate versions displayed in Fig. 8. 
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^{V2) = 



l{vi) = 



mia£{v) = 1 - u 



Fig. 8: (a) The well-labeled map with two faces and a marked vertex 
coding a triply-pointed quadrangulation with prescribed values of di3, 
and /i23, in the case /123 = (iis < ^23- The marked vertex vi is incident to 
the face /a and is connected to the frontier C2 between /2 and fs by a bridge 
b (whose edges are adjacent to /s only). The label of vi is and the minimal 
label for vertices incident to /2 (respectively /a) is 1 — t = 1 — d23 + ^123/2 
(respectively 1 — u = 1 — /123/2). The minimal label on the cycle C2 is 0, 
as well as that on the bridge b. (b) The well-labeled tree with two marked 
vertices coding a triply-pointed quadrangulation with prescribed values of 
di3, d23 and Z123, in the case /123 = dis = d23. The marked vertices vi and 
f 2 are connected by a branch b and have label 0. The global minimal label 
is 1 — w = 1 — /123/2, while the minimal label on the branch b is 0. 



2.2. Generating functions 

Known generating functions 

We can now readily relate the generating functions of the various well-labeled maps 
above to those introduced in Ref. |T^. As usual, we attach a weight g to each edge of 
a well-labeled map, which amounts to a weight g per face of the quadrangulation. The 
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i{v) = £ 




£{v) > 1 - s 




i{v) = £ 



^}) y y-) 




e{v) >i-t 




Fig. 9: A schematic picture of the known generating functions Rg, Xg^t, 
Xe-s,t and Ys^t,u (see text). 

first generating function is that of weU-labeled trees planted at a corner with label i > 
and whose labels are all larger than or equal to 1 (see Fig. 9). It reads 



Rf = R 



[i + ll [i + 2i 



where 



and where 



X 



(2.10) 



(2.11) 



R 



X 



1- y/T^n^ 

6g 

l-2Ag- Vl - 12^^ + v^7VT%7^^W^ 

2(%T7r^W^^T) 



(2.12) 
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Note that Ri = 1 + 0{g) for all £ > 1, with a conventional weight 1 for the tree reduced 
to a single vertex. The generating function of well-labeled trees planted at a corner with 
label i > and whose labels are all larger than or equal to 1 — s, for some s > 0, is then 
simply given by Ri+s, as obtained by a simple shift of all labels by s. 

The second generating function is that of well-labeled trees with two distinct marked 
vertices having label 0, connected by a branch with no n- negative labels only, and such 
that the trees attached to one side of the branch have labels larger than or equal to 
1 — s and those attached to the other side have labels larger than or equal to 1 — t, with 
s > and t > (see Fig. 9). By convention, the trees attached to the marked vertices 
are assumed to be on opposite sides, so that the result is symmetric in s and t. This 



generating function reads [[15 



m — 1 

Xs,t = X X n gRek+^Rik+t 

m>0 M = (o=to.ei,-..,em.=o) k=0 /r, i o\ 

s.t. fi>0, Ki + l-^il<l, i = 0,...,m-l (Z.io) 

[3]. [^ + iL [t + Mx [s + t + 3i 
[iL[« + 3L [^ + 3L [s + t + il- 

Note that Xs^t = 1 + C>{g), with a conventional weight 1 for the tree reduced to a single 
vertex, which is added for convenience to the family of trees enumerated by Xs^t- 

We may instead consider well-labeled trees with two marked vertices, one with label 
£ > 0, the other with label 0, with strictly positive labels on the branch inbetween and 
such that the trees attached to one side have labels larger than or equal to 1 — s and 
those attached to the other side have labels larger than or equal to 1 — t, with s > 
and t > (see Fig. 9). The tree attached to extremity with label £ is assumed to have 
labels larger than or equal to 1 — s and that attached to the extremity with label is 
assumed to have labels larger than or equal to 1 — t. The resulting generating function 



reads \15 



m — l 



[s + ll[s + 2l[t]Jt + 3l[2£ + s + t + 3] 



m>£ j^=(i=eo,ei,....ern=o) fc=l ,^ ^ ^, 

s.t. ii>0, Ki + l-^il<l, i = 0,...,m-l (2.14) 



[s + t + 3l[£ + s + ll[£ + s + 2l[£ + tl[£ + t + 3l ' 

This last formula extends to £ = where it yields Xo-s,t = 1, corresponding again to a 
conventional weight 1 for the tree reduced to a single vertex. 

The final generating function counts well-labeled trees with three marked vertices, 
say Wi, W2, Ws, and with the following constraints (see Fig. 9). On the tree, the marked 
vertices are connected by three branches joining at a central vertex. We impose that the 
branches leading respectively to wi, W2 and ws appear clockwise around this central 
vertex. We also impose that all labels on these branches be strictly positive, except 
for wi, W2 and w^, which have label 0. We further impose that trees attached to the 
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branch from wi to W2 on the side opposite to have labels larger than or equal to 
1 — s. Similarly, we impose that trees attached to the branch from W2 to (respectively 
from ws to Wi) on the side opposite to wi (respectively W2) have labels larger than or 
equal to 1 — t (respectively 1 — u). By convention, the labels on the tree attached to wi 
(respectively W2 and ws) are assumed to be larger than or equal to 1 — s (respectively 
1 — t and 1 — u). The corresponding generating function reads : 



,t,u — ^ -^£;s,t-^£;t,u-^ t;u,s 

^=0 (2.15) 

Again, we have Ys^t.u = 1 + ^{q)-, with a conventional weight 1 for the tree reduced 
to a single vertex, which is added for convenience to the family of trees enumerated by 

^ s,t,u- 

Application to minimal separating loops via the Schaeffer bijection 

In this approach, we have to enumerate trees of the type displayed in Fig. 5. It is 
convenient to first relax the condition on the global minimum, demanding only that it 
be larger than or equal to 1 — w. We can then decompose the tree by cutting it at the 
first and last occurrence of the label on the branch from vi to V2, resulting in three 
trees counted respectively by Xs-u,m Xu,u and Xt-u,u (see Fig. 10 for an illustration). 
The corresponding generating function therefore reads: 

^s+t [3]. [u]l [u + 1]^ [u + 2]l [2s + 2« + 3], [2t + 2u + 3], (2-16) 



X 



[1]. [2« + 1]. [2^ + 3], nto [s + u + kUt + u + k]^ 



To restore the condition that the global minimal label be exactly 1 — w, we simply have 
to consider A^jii/'ioop(s, t, "u) where is the finite difference operator: 

^^f[u) = f{u)~f{u-l) . (2.17) 

To conclude, the generating function for triply-pointed quadrangulations with prescribed 
values of (ii3, ^23 and /123 is given by 

G'loop('^135 '^23; ^123) = ^uH\oop{s, t, u) 

(2.18) 

with s = (ii3 - /123/2 , t = d23- /123/2 , u = /123/2 . 



Application to minimal separating loops via the Miermont bijection 

In this approach, we simply have to enumerate maps of the type displayed in Figs. 7 
and 8. Again, we relax the conditions on the minimal label within each face, namely 
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first lastO 
I 




Fig. 10: The cutting of a well-labeled tree of the type of Fig. 5 (with a 
relaxed constraint on the global minimal label) at the first and last label 
encountered along the branch from vi to V2- This results into three pieces, 
enumerated by Xs-u,u, Xu,u and Xt-u,u respectively. 

we demand only that it be larger than or equal tol — s, 1— tori — m respectively. 
In the generic case of Fig. 7, we can now decompose the map by cutting it at the first 
and last occurrence of the label on the cycle ci, starting from the endpoint of the 
bridge 6, at the first and last occurrence of the label on the cycle C2, starting from the 
other endpoint of the bridge 6, and finally at the first and last occurrence of the label 
on the bridge b itself (see Fig. 11 for an illustration). This results in general into five 
trees counted respectively by Xs^u, ^s,m,u, Xu,u, Yt,u,u and Xt^u- The corresponding 
generating function therefore reads 



-^l00p('S7 ^5 ^) Xg^^ Ys^u,U Yt,U,U Xf^^ 

[3]. [s + 1]. [t + 1], [u + 1]^ + 2^ + 3], [t + 2u + 3], 

[1]' + ^ + 1]. + « + 3] + + 1] + + 3] + 1] + 3], 

(2.19) 

Note that this formula incorporates the cases where some of the cutting points above 
coincide as we added in Xs^t and Ys,t,u the weight 1 of the tree reduced to a single 
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Fig. 11: The cutting of a well-labeled map of the type of Fig. 7 (with 
relaxed constraints on labels inside each face) at the first and last label 
encountered along the cycles ci, C2 and the branch b (see text). This 
results into five pieces, enumerated by Xg^m Ys,u,u, -^u,it5 Yt,u,u and Xt^u 
respectively. 



vertex. It also naturally incorporates the degenerate cases of Fig. 8: for instance, the 
situation of Fig. 8-(a) is properly taken into account by having the two leftmost trees in 
the decomposition of Fig. 11 reduced to single vertices, while the situation of Fig. 8-(b) 
is properly taken into account by having the two leftmost and the two rightmost trees 
in the decomposition of Fig. 11 reduced to single vertices. Again, we can restore the 
constraint that the minimal label within each face be equal tol — s, 1 — tori — w 
respectively by considering A^A^ A^jFioop(s, t, u). We deduce the alternative formula: 

G'loop('il3, <^23; ^123) = AgAtAuFioopis, t, u) 

(2.20) 

with s = di3 - /123/2 , t = d23- ^123/2 , u = /123/2 . 
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Note that the two expressions ( 2.18 ) and (|2.20 ) are consistent as we have the 
identity: 

As + u + l]Js + u + 3l Yll^^[s + u + kl 

which can be checked directly from the definition ( J2.11D . 

A simpler generating function is that of triply-pointed quadr angulations with a 
prescribed value of Z123 only. The corresponding generating function G'ioop(^i23) is ob- 
tained by summing Giooplc^is, (^23; ^123) over all the allowed values of dis and for a 
fixed /123. This amounts to a summation over all non-negative values of s and t, which 
is easily performed upon using the expression (|2.20|) by noting that, with the above 
expression ( |2.19| ), the quantities Fioop(— 1, t, u) and -Fioop(s, —1, u) vanish identically, so 
that: 

Gloop(^123) = A„Fioop(0O, 00, u) 

[1]' 12«+1L (2" + 31. 

2.3. Continuum limit 

The scaling limit is obtained by letting g approach its critical value 1/12 and 
considering large values of dis, (i23 and /123 with the following scaling: 

^ = — (1- Ae) 

12^ ' (2.23) 

1 n -1/4 J 7-) -1/4 7 T -1/4 

"13 — -t^ise 5 "23 — -'-^23e 5 n23 — J^123^ 5 

and e — 0. The quantity A may be interpreted as a "cosmological constant". In this 
limit, we have: 

G'ioop(^i23) ~ e"-^/"* 2^ioop(-^123;a) 

, ^ ,r , 1 « 3 sinh^(a[/) 3 sinh(aLi23/2) (2.24) 

where ^ioop(i^i23; a) = 7; Ou — . n, ~ 



2 a^sinh^(2at/) 



c/=Li23/2 4acosh2(aLi23/2) 



Here and throughout the paper, we use the notation: 



a 



V^AV4 _ (2.25) 



Note the factor 1/2 in the definition of ^loop, which is introduced to compensate the 
fact that, at the discrete level, Z123 can take only even integer values. More generally, 
we have 

G'ioop('^i2, c^i3, ^123) ~ e^/^ 2 ^ioop(-Di2, -D13, -^123; a) , 

Fioop(s, t, u) ~ ^ioop(^, T, U; a) , (2.26) 

-f^ioop(s, t, u) ~ 7iioop(5', T,U;a) , 
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where 

^ , ^ ^ , 3 sinh(a>S) sinh(aT) smh^(at/) smh(a(^ + 2U)) smh(a(T + 2U)) 



a" 



( smh(a(^ + U)) sinh(a(T + U)) smh(2at/))' 



^ sinh«(«^) sinh(2a(^ + U)) sinh(2a(T + [/)) 
/tioop(,<^, J , ty ; Q!j = ^ ^ ^ . — . 4 



sinh^(2at/) sinh^(a(^ + f/)) smh^(a(T + U)) 



(2.27) 



and where 



1 1 

^ioop(^i2, ^23, ^123; a) = ;j dsdrduTioopiS, T,U;a) = - duH\oop{S, T,U;a) 

^ ^ (2.28) 

with ^ = 1^13 - L123/2 , T = D22,- L123/2 , t/ = L123/2 . 

Again the two expressions above for ^ioop(-Di2, -D23, -^^123; ol) are consistent as we have 
the identity 

/ 1 sinh(a>S) sinh(a(^ + 2t/)) \ _ ^mS^{aU) sinh(2a(>5 + t/)) 

sinh2(a(5 + t/)) ; sinh4(a(5 + t/)) ' ^' ^ 

which is the continuous counterpart of ( |2.21|) . 
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L 



123 



Fig. 12: Plot of the probabihty density pioop(-^i23) 



The above continuous formulas can be used to capture the statistical properties of 
triply-pointed quadrangulations with fixed size, i.e. with a fixed number n of faces, in 
the limit n — 00. Indeed, fixing n amounts to extracting the term of the various 
discrete generating functions at hand. This can be done by a contour integral in g which, 
at large n, translates via a saddle point estimate into an integral over a real variable ^. 
More precisely, considering for instance the generating function G'ioop(/i23), we write 



G'loop(^123)|g" — /^^G'loop(^ 

J 9 



123 



) 



(2.30) 
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and we perform the change of variables 



^^Uy-^n) ' ^123 = ^123^^/' ■ (2.31) 

At large n, the contour integral becomes at dominant order an integral over real values 
of ^ and we can use the continuous formulas above with e = 1/n and A = — After a 
proper normalization by the number of triply-pointed quadrangulations with fixed size 
n, we obtain the probability density pioop(-^i23) for the rescaled length L123: 

Pioop(i^i23) = ^ / ^e-^' <3ioop{Li2z] V-3ie/2) . (2.32) 

The quantity pioop(-^i23) (^-^123 is the infinitesimal probability that the (rescaled) min- 
imal length for loops having origin ^3 and separating vi from V2 lies in the range 
[L123, Li23 + (iLi23] in the ensemble of triply-pointed quadrangulations with fixed size n, 
in the limit n 00. This probability density is plotted in Fig. 12 and has the following 
limiting behaviors: 



3 

Pioop(i^i23) ~ Y^^i23 when L123 , 

1 2/3 _(' 3 "l^/^ r ''/^ 

Pioop(^i23) ~ ^^76^123 ^ """^ when L123 00 



(2.33) 



The associated average value of L123 reads 

(L123) = ^ (£>)= 2.36198 ■■ ■ with (D) = 2^T = 1.77U8 ■ ■ ■ (2.34) 

Here and throughout the paper, we decide to express average distances in units of the 
average distance (D) between two uniformly chosen vertices in a large quadrangulation, 
whose value given above was computed in Refs. |T^-|2^]. 

Similarly, the joint probability density for Dis = di^/'n}/'^, D23 = d^zjrt}!'^ and 
L123 reads: 

Pioop(I^i2, 1^13, ^123) = ^ / dii e-^' ^?ioop(I^i2, 1^13, i^i23; ^-31^/2) , (2.35) 

while the conditional probability density for D13 and -D23, given the value of L123, 
simply reads: 

Ploop(I^12, I^13|Ll23) = ^^-P(^^3;f^3-^^^^) . (2.36) 

Ploop(,-t^l23j 

This conditional probability density is represented in Fig. 13 for decreasing values of 
L123 (namely L123 = 2.0, L123 = 1.6 and L123 = 1.0), and in Fig. 14 for increasing 
values of L123 (namely L123 = 2.0, L123 = 3.0 and L123 = 4.0). For large enough L123, 
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Fig. 13: Plots of the conditional probability density pioop(-Di3, -D23|-^i23) 
for Li23 = 2.0, Li23 = 1.6 and L123 = 1.0, from top to bottom. For each 
plot on the left, we display its associated contour plot on the right. 



this joint probability density is maximal for equal values of Di^ and -D235 i-e- when the 
two vertices vi and V2 are equally distant from vs- On the contrary, for small enough 
-^123, we observe a symmetry breaking phenomenon with a probability density being 
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maximal when one of the two vertices vi or V2 hes closer from f 3 than the other. 
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Fig. 15: (a) Plot of the conditional probability density pioop(-Di3, -D23|-^i23) 
for a small value of -L123, here L123 = 0.02. This density is concentrated 
in two regions corresponding to either Di^ or D23 being of order L123. 
A zoom on the first region is obtained by considering the same plot (b) 
with a rescaled abscissa u = 2D13/L123, or the corresponding contour plot 
(c). As apparent by taking longitudinal and transverse cut views along the 
thick lines in (c), the probability density factorizes in this region into the 
product of the density i/^iu) (red curve in (d)) and the two-point function 
p(L>23) (green curve in (e)). 



This phenomenon increases for smaller L123 and, when L123 0, we find that 

ploop(I^13, I^23|i^l23) ~ P{D^^) X -—'iP + p(I^23) X -— t/^ — ^ (2.37) 

-t^l23 V -^123 / -t^l23 V -^123 / 

with a scaling function 

i.{u) = I ^ (2.38) 

normalized to 1/2 when uj varies from 1 to oo, and where p{D) is the so called canonical 
two-point function, which is the probability density for the distance D between two 
vertices picked uniformly at random in a large quadrangulation. This canonical two- 
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point function is given by a formula similar to (|2.32 ): 



p{D) = ^ / d^^ e-^" G{D- 7=31^72) 

'^"^ ^--^ (2.39) 



with g{D]a)=Aa''- 



3 cosh(aL>) 



sinh (tt-D) 

The particular form ( |2.37| ) expresses that, when L123 becomes small, one of two vertices 
Vi or f2, say vi necessarily lies in the vicinity of ^3, with a distance D13 of the order of 
L123 and governed by the density ( ]2.38| ) for u = 2D13/L123, while the other vertex lies 



at an arbitrary distance from the two others, with a probability density given simply 
by the two-point function of quadrangulations, as expected. This behavior is depicted 
in Fig. 15, for L123 = 0.02. This result corroborates the known property of quadran- 
gulations of large size n: small loops of length negligible with respect to n^/^ in the 
quadrangulation necessarily separate it into a macroscopic domain containing most of 
the area of the quadrangulation, and a small part of negligible size with respect to n. 



loop 






Fig. 16: (a) Plot of the conditional probability density pioop(-Di3, -D23|-^i23) 
for a large value of L123, here L123 = 20.0. The same plot (b) and its con- 
tour (c) in the rescaled variables n = (-D13 — Li23/2)(9Li23/2)^/'^ and 
ly = (-D23 — -^i23/2)(9Li23/2)^/'^. In these variables, the conditional prob- 
ability density tends to a limiting distribution $(^, z^), as shown in (d). 
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In the other hmit, i.e. when L123 becomes large, we find the hmiting behavior: 



Ploop(I^13, ^231^123) ~ ( ^ ) i^) 



or \ 



-^123 \ /9Li23\^^'^ -^123 A /9Li23^^^'^ 



(2.40) 

with n= iDis — j ( j ^ u = ID23 

with a scaUng function 

$(^,zy) =e-('^+^) (2-e-'^-e-^) (2.41) 

properly normalized to 1 when fu, and u vary from to 00. At large L123, both distances 

Dis and D23 are therefore necessarily of order L123/2, with differences Dis — L123/2 

]^ 73 

and L>23 — L123/2 of order L^23 ? governed by the joint probability density ( |2.41| ). This 
behavior is depicted in Fig. 16 for L123 = 20.0. 



3. Confluence 

3.1. Confluence of geodesies 




Fig. 17: A schematic picture of the phenomenon of confiuence of geodesies. 
For generic points vi , f 2 and f 3 and in the scaling limit of large quadran- 
gulations, the geodesic from vi to vs and that from V2 to ^3 (represented 
as thick blue lines) are unique and have a common part of macroscopic 
length 6. 
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In this section, we explain how we can use the quantities computed in previous 
section, or shght generahzations of them, to study the phenomenon of confluence of 
geodesies in the scahng hmit of large quadr angulations. It was shown by Le Gall HIT 



and Miermont [T^] that two typical points in a large random quadrangulation are joined 
by a unique "macroscopic" geodesic path. By this, it is meant that, although there is a 
large (extensive in the length) number of geodesic paths between two points at a discrete 
level, all these geodesies remain within a distance negligible with respect to n^/^, which 
is the scale at which points can be distinguished in the scaling limit. Moreover, given 
three typical vertices fi, V2 and fs, the unique macroscopic geodesic from vi to V3 
and the unique macroscopic geodesic from V2 to vs merge before reaching V3, i.e. have 
a macroscopic common part (see Fig. 17). This is the phenomenon of confluence of 
geodesies [jl^ which raises interesting problems, such as that of the distribution of the 
length 6 of this common part. 



Approach via the Schaeffer bijection 



(a) 



1{V2) = t 





--■ l-u'" 



yt-2 



(b) 



Fig. 18: In the well-labeled tree of Fig. 5, we distinguish (a) the minimal 
label 1 — u' on trees attached to one side of the branch from vi to f 2 and 
the minimal label 1 — u" on trees attached to the other side of the branch, 
with u = max(w', u"). The quantity \u' — u"\ measures the length of the 
common part of the leftmost geodesies from vi and f 2 to the added vertex 
V3. As apparent in (b), here in the case u' > u" , these leftmost geodesies 
are made of two distinct chains of successors of respective lengths s — u" 
(green long-dashed arrows) and t — u" (red short-dashed arrows) , followed 
by a common chain of successors of length u' — u" (magenta solid arrows). 



At a discrete level, this length can be estimated by a particular choice of geodesies 
defined as follows: we start again with a triply-pointed quadrangulation with marked 
vertices vi, V2 and ^3 and consider the associated well-labeled tree obtained from the 
Schaeffer bijection, taking ^3 has the origin. This tree has two marked vertices vi and 
f2, and upon shifting the labels so that the minimal label on the branch between vi 
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and V2 is 0, it is of the type displayed in Fig. 5 for some s, t and u. We can now 
consider the leftmost geodesic from vi to f 3 formed by the chain of successors from 
the corner incident to Vi and lying immediately on the right of the branch oriented 
from vi to V2- Similarly, we consider the leftmost geodesic from V2 to obtained as 
the chain of successors from the corner incident to f 2 and lying immediately on the 
left of the branch (oriented again from vi to ^2). These two geodesies will merge at a 
point which we characterize as follows (see Fig. 18 for an illustration): let us call 1 — u' 
(respectively 1 — u") the minimal label on trees attached to the left (respectively right) 
side of the branch oriented from vi to f 2 (with the convention that the tree attached 
to Vi lies on the left side of the branch, and that attached to V2 on the right side), 
with u = max('u',u"). Then the two chosen geodesies have a common part of length 
\u' — u"\. Indeed, assuming without loss of generality that u = u' > u" , all the s + u 
successors of the corner chosen at vi lie on the left of the branch until f 3 (with label 
—u) is reached. On the other hand, among the t + u successors of the corner chosen at 
f 2, the first t + u" — 1 successors are found on the right of the branch but the [t + u")-th. 
successor, having label —u", is on the left of the branch and coincides with the {s + u")- 
th successor of the corner chosen at vi. From that point, all remaining successors form 
a common part of length u' — u" (see Fig. 18). To conclude, there is a correspondence 
between, on the one hand, well-labeled trees with fixed values of s, t, u' and u" as 
defined above and, on the other hand, triply-pointed quadrangulations with prescribed 
values (ii3 = s -f max(M', u"), ^23 = t + max('u', u"), /123 = 2max(w', u") and such that 
the leftmost geodesies from vi to f 3 and from V2 to f 3 have a common part of length 
\u' — u"\. Note that the sign of u' — u" simply accounts for the relative position of the 
geodesies: when u' > u" (respectively u" > u'), the geodesic from vi to ^3 merges on 
the right (respectively on the left) of the geodesic from V2 to ^3. 

We now wish to enumerate the above trees. By an immediate generalization of 
Eq. (|2.16| ), such trees have generating function: 

A^,/ A„//ifioop(s, t, u , u") where ifioop(s, t, u', u") = Xs.u',u" Xu',u" Xt-u",u' ■ (3.1) 
In the scaling limit, this generating function becomes: 

du'dvUioo^iS, T, U', U"; a) where 
nioop{S,T, U',U";a) = 

g sm\i^{aU') sinh^(a;7") smh{a{2S+U' + U")) smh{a{2T+U' + U")) 

{smh{a{U' + U")) sinh(a(^+C70) smh{a{S+U")) smh{a{T+U')) smh{a{T+U")))^ 

(3.2) 

and we expect that any other choice for the geodesies at the discrete level would lead to 
the same continuous expression. This formula holds in the grand-canonical formalism 
and can be transformed via an integral of the type ( p.32| ) into the canonical normalized 
joint probability density for -D13, -D23, -^123 and the (rescaled) length S = \U' — U"\ for 
the common part of the geodesies. 
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(b) 

Fig. 19: In the well-labeled map of Fig. 7, we mark the last occurrence 
of a label on each (counterclockwise oriented) cycle ci and C2 and call 
1 — u' the minimal label on trees attached to the part of the frontier of 
the external face made of: (i) the left side of the branch b (oriented from 
ci to C2), (ii) the external side of the cycle ci before reaching the marked 
label 0, and (iii) the external side of the cycle C2 after passing the marked 
label 0. We also call 1 — u" the minimal label on trees attached to the 
complementary part of the frontier, with u = m.ax{u',u"). The quantity 
\u'—u"\ measures the length of the common part of two particular geodesies 
leading from vi and V2 to ^3, as apparent in (b), here in the case u' > u" . 

Approach via the Miermont bijection 

As in section 2, a useful alternative expression for the above function may be 
obtained by use of the Miermont bijection for triply-pointed quadrangulations, leading, 
for the special choice ( p.8| ) of delays, to well-labeled maps of the type displayed in Fig. 7 
(or of its degenerate versions) for some s, t, and u. A particular geodesic path from vi 
to f 3 is obtained by picking say, the last label on the (counterclockwise oriented) cycle 
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Ci , looking at the two corners at that vertex lying immediately on the right of the cycle 
when we follow the cycle in both directions, and considering the chains of successors of 
these two corners. The concatenation of these chains forms the desired geodesic path. 
A similar geodesic path can be considered from V2 to fs, passing via the last label on 
the (counterclockwise oriented) cycle C2- Let us now call 1 — u' the minimal label on 
trees attached to the left side of the branch b (oriented from vi to ^2), to the external 
side of the cycle ci before the last occurrence of a label on this cycle, and to the 
external side of the cycle C2 after the last occurrence of a label on this cycle (see 
Fig. 19 for an illustration). We also call 1 — u" the minimal label on trees attached to 
the complementary part of the frontier of the face /s, with u = max(w', u"). Then by 
arguments similar to the discussion above, the two particular geodesies have a common 
part of length \u' — u"\. We now have a correspondence between, on the one hand, 
well-labeled maps with fixed values of s, t, u' and u" as defined above and, on the other 
hand, triply-pointed quadr angulations with prescribed values di^ = s -|- max(u', u"), 
<^23 = t + max(M', w"), /i23 = 2max(u', w") and such that the two particular geodesies 
considered above from vi to ^3 and from V2 to have a common part of length \u' — u"\. 
By an immediate generalization of Eq. (|2.19[) , such maps are enumerated by: 



/\s^t^u'^u"F\oop{s,t,u ,u ) where 

(3.3) 

-^loop('S, t,U ,U ) = Xs,u' Ys,u",u' Xu',u" Yt,u',u" Xt^u" 



Note that this generating function is different from that given by (|3.1| ) as our particular 
choice of geodesies differs in the Schaeffer and in the Miermont bijection approach. In 
the scaling limit however, we expect to recover the same expression ( |3.2| ) due to the 
unicity of geodesies at a macroscopic level. Indeed, the expression translates into: 



dsdrdu'du" J^ioopiS, T, U', U"\ a) where 
J'ioop(5', T, U' , U"; a) = 

3 sinh(a^) sinh(aT) sinh2(at/') smh^iaU") smh{a{S+U' + U")) smh{a{T+U' + U")) 

smh{a{S+U')) smh{a{S+U") smh{a{T+U')) smh{a{T+U")) smh^{a{U' + U")) 

(3.4) 

which precisely matches the continuous expression (|3.2| ), namely: 



dsdTdu'du"J'ioop{S, T, U', U"; a) = du'dwHioopiS, T, U', U"; a) (3.5) 
as a consequence of the identity 

„ fl sinh(a^)sinh(a(^+t/' + t/")) \ sinh(at/') sinh(aC/") sinh(a(2^+C/' + t/")) 



a sinh{a{S+U')) siuh{a{S+U")) J sinh^ {a{S+U')) sinh^ {a{S+U")) 

(3.6) 
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Marginal law for S 

It is now a simple exercise to obtain, in this scaling limit, the marginal law for 
5. We simply have to integrate over all positive values of S, T, U' and U" with the 
constraint that \U' — U"\ = 6. This is done more easily in the grand-canonical formalism 
first and by use of the expression (|3.4|), namely: 



l*OC i*oo i*oo i*oo 

/ dS dT dU' / dU"6{\U' -U"\- 6) dsdTdu'du"J^ioop{S, T, U', U' 
Jo Jo Jo Jo 

/oo POD 
dU' / dU"6{\U' -U"\- 6) du'du" J^ioopioo, oo, U\ U") 

= rciU' hu'm' - U"\ - 5) Su.Ou.. f4^'"'"(aC/-)sinh'(aC/")X 
Jt, Jo \ I u u smh2(a(f/' + f/")) / 

= / dU' dU"5{\U' -U"\-5) 18 
Jo Jo 



36 / dU 



sinh^(a(t/' + U")) 
sinh2(aC/) sinlv^ {a{U - 5)) 



sinh^(a(2f/ - 5)) 



2a 



(3.7) 



a(5) 




Fig. 20: Plot of the probability density a(5) for the length 5 of the common 
part of the two geodesies from vi and V2 to in the scaling limit of large 
triply-pointed quadrangulations. 



As before, we can transform this result into the probability density for the (rescaled) 
variable 5 in the canonical ensemble of triply-pointed quadrangulations of large fixed 
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size. This probability density reads: 



a{5) 



i\/n 



2a 



1 1 

4' 2 



a=^-3i5/2 



' 64 ^ 
3 31 

4' 2 J '""el 

5 71 



4' 4 



64 



where 



00 I, 



1 



fc-1 



k 



with = n^^+' 



(3.8) 



(3.9) 



fc=0 " i=0 

This probabihty density is plotted in Fig. 20. We have in particular 

(S) = 1 (D) = 0.590494 ■ ■ ■ 

i.e. the common part represents on average one third of the length of a geodesic. 
3.2. Confluence of minimal separating loops 



(3.10) 




Fig. 21: A schematic picture of the phenomenon of confluence of mini- 
mal separating loops. For generic points vi, V2 and V3 and in the scaling 
limit of large quadr angulations, the minimal loop originating from f 3 and 
separating vi from V2 (represented as a thick blue line) is unique and is 
made of a common part of macroscopic length 5ioop and an open part of 
macroscopic length L^23- 
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The minimal separating loops themselves also exhibit a phenomenon of confluence. 
Indeed, a minimal separating loop is made of two geodesies of same length emanating 
from a particular vertex v (with minimal label on the branch from vi to f 2) and reaching 
Vs. In the scaling limit, we expect that the macroscopic minimal separating loop is 
unique and moreover, its two constituent geodesies have a common part of macroscopic 
length 5ioop (see Fig. 21 for an illustration). Note that, although vi, V2 and ^3 are 
generic points, f is a non-typical point as it can be connected to by two distinct 
macroscopic geodesies which are not confluent at v. We shall call the complementary 
part the open part of the loop, with length L']^23 = -^123 — 2(5ioop- 




(b) 

Fig. 22: In the well-labeled tree of Fig. 5, we mark the flrst label on the 
branch from vi to V2 and call 1 — u' the minimal label on trees attached 
to the part of the branch lying from vi to the marked label 0. We also 
call 1 — u" the minimal label on trees attached to the complementary part 
of the branch, with u = max(u',zi"). The quantity \u' — u"\ measures 
the length of the common part of a particular minimal loop originating 
from f3 and separating vi from ^2, as apparent in (b), here in the case 
u' > u" . The length of the open part of the minimal separating loop is 
2min(u', w") = 2u". 

The statistics for 5ioop and L'123 can be computed along the same lines as in the 
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previous Section. In the Schaeffer approach, on the branch from vi to f 2 in the well- 
labeled tree, we now consider the vertex v with minimal label closest to vi. Calling 
1 — u' the minimal label for trees attached the the part of the branch from vi to v, and 
1 — u" the minimal label for trees attached to the complementary part, the quantity 
\u' — u"\ measures the length of the desired common part for a particular minimal loop 
formed by two chains of successors starting from v (see Fig. 22 for an illustration). As 
for the length of the open part of the loop, it is simply measured by 2min(w', u"). The 
generating function for the objects above is immediately given by 

Au'Au"Hioop{s, t, u\ u") where -ffioop(s, t, u', u") = Xs;u',u' Xu",u" Xt-u" ,u" ■ (3.11) 

In the scaling limit, this generating function becomes: 

du'du"H\oop{S, T, U', U"; a) where 

^ ^^rj^rr'TT" ^ , sinh"(«^0 sinh"(«^-0 sinh(a(2(^+^0) sinh(a(2(r+[/-0) 
^ sinh(2aC/0sinh(2af/")sinh^(a(5+f/'))sinh^(a(T+t/") 

(3.12) 

which yields the joint law for L'123 = 2min(t/', U"), Dis = S + max(t/', U"), D23 = 
T + max(t/', U") and Sioop = \U' -U"\. Note that the sign of U' - U" indicates which 
domain delimited by the open part contains the common part of the loop. As apparent 
in Fig. 22, this common part lies in the domain containing vi when U' > U". 

An alternative expression is found through the Miermont approach where we con- 
sider well-labeled maps of the type displayed in Fig. 23 using a particular minimal 
separating loop passing through the vertex with minimal label on the branch b closest 
to the cycle ci. We find a generating function 



AsAtAu'Au"Fioop{s, t, u\ u") where 

-^loop('S, u' , u") = Xgy Ysy^u' Xu",u" Yt,u",u" Xt^ 



(3.13) 



whose scaling limit 



dsdTdu'du"J^\oop{S, T, U', U"; a) where 
-^ioop(<S', T, U',U";a) = 

3 sinh(ag) sinh(aT) smh^jaU') sinh^iaU") sinh(a(^+ 2U')) sinh(a(r+ 2U")) 
sinh^{a{S+U')) sinh^ {a{T+U")) sinh(2at/0 sinh(2a[/") 

(3.14) 

matches the expression ( p.l2| ) above. This matching is again a direct consequence of 



the identity (|2:29|) . 



We can integrate over S and T and obtain the marginal law for U' and U" 

Q sinh'{aU')sinh^{aU") _ 3 1 

a2 ^' ^" sinh(2ctC/0 sinh(2aC/") 4. cosh^{aU') cosh^{aU") ' 
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Fig. 23: In the well-labeled map of Fig. 7, we mark the first occurrence 
of a label on the branch b oriented from the cycle ci to the cycle C2 and 
call 1 — u' the minimal label on trees attached to the part of the frontier 
of the external face made of: (i) the part of the branch h lying between ci 
and the marked label and (ii) the external side of the cycle ci. We call 
1 — u" the minimal label on trees attached to the complementary part of 
the frontier, with u = max(w',w"). The quantity \u' — u"\ measures the 
length of the common part of a particular minimal loop originating for f 3 
and separating vi from V2, as apparent in (b), here in the case u' > u" 
(corresponding to having the common part in the domain containing Vi). 
The length of the open part is 2min(u',u"). 



This can be translated into the marginal joint law for 5ioop and ^^'123- For triply-pointed 
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Fig. 24: Plot of the joint probability density rioop('^ioop, -^123) ^^e 
length 5ioop of the common part and the length L']^23 of the open part of 
the minimal loop originating from and separating vi from V2 in the 
scaling limit of large triply-pointed quadrangulations. 

quadrangulations of fixed large size n, we find the joint probability density: 



''"loopl^loopi -^123) 



3/4 



cosh^ fttf -^^+5ioop) ) cosh^ 



This probability density is plotted in Fig. 24. 



V-3i?/2 

(3.16) 




loop 



Fig. 25: Plot of the marginal probability density rioop(5ioop) for the length 
of the common part of the minimal loop originating from f 3 and separating 
vi from V2 in the scaling limit of large triply-pointed quadrangulations. 
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Upon integrating rioop('Jioop7 -^123) o"^^^ -^'1235 S^^ marginal density distribu- 
tion for the length 5ioop only, namely: 



noop(5: 



loop J 



U/TT 



2a sinh^(a(5ioop) cosh(a;5ioop) 
X <; 2cosh^(Q!5ioop) (a5ioop-log(cosh(a5ioop))) 

-sinh(a5ioop) (cosh(a5ioop) +e""'^'°°p) 



(3.17) 



a=^-3ie/2 

This probability density is plotted in Fig. 25. We have in particular 

2 



(5loop) 



(2 -log 4) (D) = 0.724779- 



(3.18) 



in terms of the average length {D) of a geodesic path. 

T (r ' ) 




^123 



Fig. 26: Plot of the marginal probability density Tioop(-^i23) length 
of the open part of the minimal loop originating from f 3 and separating vi 
from f 2 in the scaling limit of large triply-pointed quadrangulations. 

On the other hand, upon integrating rioop(5ioop, L'j^23) over 5ioop, we get the 
marginal density distribution for the length L'123 only, namely: 



Tloop(-t'l23) 



U/TT 



1 



^'-^ cosh^ ( a^^Y^ 



(3.19) 



a=^-3i?/2 

This probability density is plotted in Fig. 26. We have in particular 

(L;23) = ^ (log4 - 1) (D) = 0.912418 ■ ■ ■ . 



(3.20) 



Note that 2(5ioop) + (^'123) = U^), in agreement with Eq. (|2.34|) . 
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3.3. Area enclosed by a minimal separating loop 

Within the above framework, we may easily address the question of the dispatch- 
ing of the area of triply-pointed quadrangulations in the two domains separated by a 
minimal separating loop. For a given separating loop, we may indeed decide to attach 
a weight per face of the quadrangulation depending on which domain it lies in. In the 
equivalent Miermont picture and for the particular minimal separating loop considered 
in Fig. 23, this amounts to assigning a weight, say gi (respectively (72) to edges lying in 
the domain containing vi (respectively V2), which results in the generating function: 

A^AtA^^/ A„//Fioop(s, t, u', u"; gi, 5(2) where 

-f^loop('S, t, u\ u"; gi, (72) = Xs^u'{gi) Ys,u' ,u'{,gi) Xu",u"{92) Yt,u",u"i92) Xt^u"{92) ■ 

(3.21) 

Here Xs,t(fi'm) and Ys^t,u{gm) denote the generating functions Xg^t and ys,t,u5 as given 
by ( p.l3|) and (|2.15| ), with g replaced by gm in ( |2.12| ), for m = 1,2. In the continuum 



limit, we set 

^7i = ^(l-Aie) , ^72 = ^(1-A2e) (3.22) 

which amounts to having a different cosmological constant in both domains. The gen- 
erating function above translates into the scaling function 

dsdrdu'dw'f'ioopiS, T, t/', t/"; ai, 02) where 

/•ioop(^, T, U', U"; ai, ^2) = ^^y{S, U', U'; ai)y{T, U", U"; ^2) ^3 23) 

1 sinh(aS') sinh(aT) sinh(at/) sinh(ct(S' + T + U)) 



y{S,T,U;a) 



3a sinh(a(5 + T)) sinh(a(T + U)) sinh(a(t/ + S)) 



and am = -\/3/2Am^- Here, y is the scaling limit of Y, while each X tends to 3 in the 
scaling limit, irrespectively of its arguments. 

Upon integrating over all possible values of 5", T, U' and U" , we get a function 

3 

^ioop(oo, 00, 00, 00; ai, 0^2) = -. • (3.24) 

4aia2 

Returning to the canonical formalism where we fix the sizes of the two domains separated 
by the minimal loop to be respectively ni and n2 (with ni + 712 = n, ui ^ 1, 712 ^ 1), 
we set e = 1/n and A^, = — so that the expression ( |3.24| ) tends to n^/^ i/{2^/^l^). 



Setting Til = rjn and n2 = (1 — ?7)n, we obtain the probability density for rj as 



n POO />00 

1^ J —00 J —00 



{i-v)e2. 1 



(3.25) 

TT 1 



r(i) ry3/4(i_^)3/4 • 
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The dispatching r] of the mass is therefore governed by a simple Beta distribution with 
parameters {1/4, 1/4}. In particular the two domains are most likely of very asymmetric 
sizes, with a probability density maximal for 77 = or 1. 

More precisely, we can naturally distinguish the two domains as exactly one of them 
contains the common part of the minimal separating loop. As mentioned above, this 
information is encoded in the sign of U' — U" . We may integrate ( |3.23| ) over S, T, U' 
and U" in the domain U' > U" , corresponding to the case where the common part lies 
in the domain containing vi. This leads to 

dU'-^— [du' tanh(ait/')] tanh(a2t^') , (3.26) 
4aia2 

which, together with the symmetric contribution from the domain U" > U' (obtained 
by exchanging ai and 02)5 adds up to (|3.24|) . We can in principle deduce from ( |3.26|) 
the (now asymmetric) law for 77 conditionally on the position of the common part. We 
have not found a compact simple form for this law but its first few moments can be 
computed. We find an average value 

{V)w>u'' = (1 - V)u">u' = ^(1 + log 4) ~ 79.543% (3.27) 

for the proportion of the total area lying in the same domain as the common part of 
the minimal separating loop. 



4. The three-point function revisited 

The three-point function of planar quadrangulations enumerates quadrangulations 
of the sphere with three marked vertices vi, V2 and vs at prescribed pairwise distances 



di2, d23 and d^i. It was computed in Ref. [|T5| and, in the scaling limit of quad- 
rangulations of fixed large size n ^ 00, translates into a universal joint probability 
p(L>i2, -D23, -D31) for the three rescaled lengths D12 = di2/n^^^., -D23 = d2zlrt}'l^ and 
-D31 = dzxjr}}!^ of the three geodesies forming the triangle {y\,V2-,V'i). As mentioned 
in the introduction, a full description of the geometry of this triangle must incorporate 
the phenomenon of confiuence. We call the lengths of the common parts respectively 8\ 
(for the two geodesies leading to f 1), ^2 (for the two geodesies leading to vq) and ^3 (for 
the two geodesies leading to fa). The remaining proper parts of the geodesies form an 
open triangle with sides of respective lengths D'y2, = -D12 — 5i — ^2, -D23 = -^23 — ^2 — ^3 
and L>32^ = D31 — ^3 — 5i (see Fig. 27 for an illustration). A natural question is that of 
determining the corresponding joint probability density p{D[2, -D23, D'^i, ^i, ^s)- 

As explained in Ref. [^, triply-pointed quadrangulations with prescribed values 
of the pairwise distances di2, d23 and dsi are in one-to-one correspondence with partic- 
ular well-labeled maps with three faces. This is again a consequence of the Miermont 
bijection with three sources fi, V2 and ^3, and with a particular choice of delays, now 
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V3 



Fig. 27: A schematic picture of the phenomenon of confluence of geodesies 
for the three geodesies Unking three generic points vi, V2 and vs in the 
scaling limit of large quadr angulations. These geodesies (represented as 
thick blue lines) have a common parts of macroscopic lengths 5i , 82 and ^3 . 
The remaining open part of the triangle has sides of macroscopic lengths 
-D^2 5 -^23 -D31, and separates the quadrangulation into two domains. 
The three common parts may lie in the same domain (as represented here) 
or in different domains, giving rise to eight possibilities for the relative 
position of the three geodesies (see Fig. 33 below). 

given by 

Ti = -S = 



d23 


- dsi 


- di2 




2 




d3i 


- di2 


- d23 




2 




di2 


- '^23 


- dsi 



r, = -t^ , (4.1) 

T3 = -u 



2 

where we use the parametrization of the pairwise distances 



di2 


= s + t , 


d23 


= t + u, 


dsi 


= u + s 



(4.2) 



As shown in Ref. [T^], the maps obtained for this choice of delays are now of the type 
displayed in Fig. 28, or degenerate versions of this generic form when one of the frontiers 
between faces or one of the faces reduces to a single vertex. By an simple decomposition 
of the map in five pieces obtained by cutting the map at the first and last occurrence of 
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first 




last 

Fig. 28: Structure of the well-labeled maps with three faces coding triply- 
pointed quadr angulations with prescribed values of the pairwise distances 
between the marked vertices: di2 = s + t, d23 — t + u and dsi = u + s. 
These maps are easily enumerated by cutting them at the first and last 
label on each frontier (big green dots). 

a label on each frontier, we immediately get the generating function for these maps: 

AsAtAuF(sA,u) where 

(4.3 

F{s,t,u) =Xs,tXt,uXu,s {Ys,t,uf . 



If we now consider, say the first label on each of the three pairwise frontiers 
between faces, the global frontier of each face can then be divided into two parts lying 
inbetween the two marked labels on this frontier. We may then distinguish the minimal 
label on trees attached to the first part of the frontier from that on trees attached to 
the second part of the frontier (see Fig. 29- (a) for an illustration). For instance, the 
minimal label 1 — s in the first face corresponds to a minimal label 1 — s' on one part and 
1 — s" on the other part with s = max(s', s"). We have similar minima 1 — t', 1 —t" and 
1 — u' , 1 — u" in the other faces. The quantities \s' — s"\, \t' — t"\ and \u' — u"\ measure 
the lengths of the pairwise common parts of three particular geodesies made of chains of 
successors of corners at the marked labels (see Fig. 29-(b)). Similarly, the quantities 
min(s', s") -f- min(t', t"), min(t', t") -f min(w,', u") and min(u', u") -(- min(s', s") are the 
lengths of the proper parts of the same three geodesies. With the above refinements. 
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(a) 



• (b) 



Fig. 29: In the map of Fig. 28, we mark the first label on each frontier. 
We then call 1 — s', 1 — s", 1 — i', 1 — t", 1 — u' and 1 — u" respectively 
the minimal label on trees attached to the six frontier sides delimited by 
these marked points as shown, with s = max(s',s"), t = max(t',t") and 
u — ma,x{u' ,u"). The quantities \s' — s"\, \t' — t"\ and \u' — u'\ measure the 
lengths of the common parts (represented by solid magenta arrows) of three 
particular geodesies obtained from the concatenation of chains of successors 
of the marked labels 0. The situation represented here corresponds to 
s" > s', t" > t' and u" > u'. 

the generating functions now reads 

A^/ A^// At/ At// A„/ A^//F(s', s", t\ t" , u' , u") where 

F(s', s", t', t" , u' , u") = Xs"^t" Xt"^u" Xu'i^s" Ys',t',u' Ys",t" ,u" ■ 

and its continuous counterpart reads 

ds'ds"dT'dT"du'du"J^{S', S", T', T", U', U"; a) where 

J^{S', S", r, T", U', U"- a) = 333^(^', T', U'- a)y{S" , T", U"; a) 

1 sinh(a5') sinh(aT) sinh(aC/) sinh(Q;(5' + T + U)) 



y{S,T,U-a) 



(4.4) 



(4.5) 



3a sinh(a(5 + T)) sinh(a(T + U) sinh(a(t/ + S)) 



which yields directly the joint law p(D^2 5 -^235 -^si? ^i? '^2, (^s) for -D^2 = min(5",5'") + 
min(r, T"), = min(T', T") + min(t/', U"), D31 = min([/', U") + min(^', S"), 5i = 
1^' - ^"1, ^2 = \T'-T"\ and d3 = \U'-U"\. Note that the sign of S' - S" (respectively 
T' — T" , U' — U") indicates in which of the domains delimited by the open part of the 
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Fig. 30: Plots of the conditional probability density 6'(5i, 52I-D12) for the 
lengths of the two common parts of a geodesic of fixed length -D12, here 
for (a) 1^12 = 1.0, (b) 2.0 and (c) 5.0. 



triangle the common part leading to vi (respectively V2-, fs) lies. Let us now discuss in 
more details a number of marginal laws inherited from p{D[2, -D23, D'^i, <5i, ^2, ^s)- 

A first marginal law is that for the lengths 61, and 62 of the three parts of 



the geodesic between vi and V2. It is obtained by first integrating (|4.5|) over U' and 
U", which yields ds'ds"dT'dT"J^iS' , S" ,T' ,T" , 00, 00), then integrating over S', S", 
T' and T" with fixed values min(S", S") = a, max(S", S") — a + 5i, min(S", 5"') = r, 
max(S", 5"') = r+52, and finally integrating over a and r with the condition a+r = -D^2- 
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We obtain the grand-canonical function 

1 



3 

-a 
2 



+ 



sinh^(aL''i2) sinh^(a(L';2 + '^i + ^2)) sinh^(a(L>;2 + '^i)) sinh^(a(L';2 + ^2)) 
X |2Q;L'i2 (2 cosh(a5i) cosh(a52) + cosh(a(2L'i2 + (^i + ^2))) 

+ 2sinh(a(5i + 52))-2sinh(a(2D;2 + 5i + '^2))-cosh(a(5i-52))sinh(2aD;2)} 

(4.6) 

from which we can get the canonical joint probability density ^2, D'vi) as before. It 
is interesting to consider this probability density conditionally on the value of the total 
length Dyi of the geodesic between v\ and V2, namely: 

^(5l,52,I^12-5i-52) 



(4.7) 



P(^12) 

where p(-D) is the canonical two-point function (|2.39| ). This conditional probability 
density is plotted in Fig. 30 for Dyi = 1.0, 2.0 and 5.0. 
At large -D12, it takes the simple form 

1 



D 
1 



X 



12 



(9I)i2)^/3x((9I)i2)'/'(I^i2 - 5i - 52)) 



where xW = 77 



1 



+ 



3 Vsinh"'(A/2) e 



3A/2 



(Acosh(A/2) -2sinh(A/2)) 



(4.^ 



In this limit, the geodesic consists mainly of two common parts linked by a small open 

— 1/3 

part whose length is of order 5 with a distribution given by the scaling function 

x(A). The position of this open part is moreover uniform along the geodesic. This 
property is illustrated in Fig. 31. 



QiD' ) 




0.5 



Fig. 32: Plot of the probability density 9{D'^2) fo^ the length of the 
proper part of the geodesic between vi and f 2 in the scaling limit of large 
triply-pointed quadrangulations. 
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e(8i,8 




(b) 



X 



(c) 



Fig. 31: The conditional probability density 9{Si, S2\Di2) for a large value 
of Di2 (here D12 — 10.0) becomes uniform in the "transverse" direction 
(corresponding to fixing the value of 61+62) and characterized by the scal- 
ing function x(A) in the "longitudinal" direction (corresponding to vary- 
ing the value of 61 + 62) , with a scaling variable A = {9Di2y^^D[2 = 
{9Di2Y'^{Di2-6i-62). 

Upon integrating ( [4.6| ) over 61 and ^2, we can get get the marginal law for D'12 
only. In the grand canonical formalism, it reads: 



IGctsinh {aD^2 



. (^1 _ (20 - 3e-2«^- + } 



(4.9) 



from which we obtain the canonical probability density 6{D'i2)- This probability density 
is plotted in Fig. 32. We have in particular 



(L>i2) = ^(i^) = 0.590494. 



(4.10) 
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i.e. the length of the open part represents on average one third of the length of a geodesic, 
in agreement with (|3.10| ). Upon integrating (|4.6|) over 1)^2 ^2 and turning to the 
canonical formalism, we can recover the marginal law a (Si) of Section 3.1. Similarly, 
upon integrating ( [4.6|) over D'12, di and 62 with a fixed value of D12 — D'12 + Si + 62, 
and upon turning to the canonical formalism, we recover the two-point function p{Di2), 
as it should. 
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Fig. 33: A schematic picture of the eight possible arrangements for the 
three common parts of the geodesies with respect to their open part. In the 
cases (a) and (h), the three common parts lie on the same side of the open 
part. In the remaining cases, two of the common parts lie on the same side 
and the third one on the other side. Note that, due to the orientation of 
the sphere, (a) and (h) (respectively (b) and (e), (c) and (f), (d) and (g)) 
can be distinguished through the cyclic order of the three sources. 

To conclude this section, let us finally discuss the global arrangement of the three 
geodesies on the sphere. As illustrated in Fig. 33, there are eight possibilities: in two 
cases ((a) and (h) in Fig. 33), the three common parts lie in the same domain, while in 
the remaining six cases ((b) to (g) in Fig. 33), two of the common parts lie in the same 
domain and the third one in the other domain. 

We may wonder what is the probability to observe a given arrangement in the 
canonical ensemble. Any of the above arrangements corresponds simply to a choice 
of sign for S' — S", T' — T" and U' — U" . To obtain, say arrangement (a), we may 
integrate (|4.5| ) with the conditions S' = max(S", S") = S, T' = max(T', T") = T, 
U' — max(t/', U") — U, leading to the grand-canonical function 

33 {dsdrduyiS, T, U; a)) y{S, T, U; a) . (4.11) 

To obtain arrangement (b), the third condition must be replaced by U" = max(C/', U") = 
U, leading to a grand-canonical function 

3^{dsdTy{S,T,U;a)) {duy{S,T,U; a)) . (4.12) 
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All the remaining arrangements follow by symmetry and their contributions add up to 
the grand-canonical three-point function 



3''dsdTdu{y{S,T,U;a)y , (4.13) 



which is the continuous limit of Eq. (^]^) . To obtain the probability of having a given 
arrangement, we simply have to integrate its individual contribution over 5", T, U, and 
divide by the integral of ( 4.13 ). Note that this ratio, obtained in the grand-canonical 



ensemble yields directly the correct canonical probability since all grand-canonical in- 
dividual contributions integrate to a numerical constant times the same function 
A simple calculation shows that each of the arrangements (a) and (h) occurs with a 
probability 1/4, while each of the arrangements (b)-(g) occurs with a probability 1/12. 

As for the dispatching of the area in the two domains, we find that, if we disregard 
the particular arrangement at hand, the probability density for the proportion 77 of 
the total area lying in one of the two domains is again given by the symmetric Beta 
distribution ( |3.25| ) with parameters {1/4,1/4}. On the other hand, if we consider a 
particular arrangement, the dispatching of the area in no longer symmetric in the two 
domains. In the case of arrangement (a) or (h), we find that, on average, ~ 94.259% of 
the total area lies in the domain containing the three common parts, while, in the case 
of the arrangement (b), (c), (d), (e), (f) or (g), an average of ~ 67.224% of the total 
area lies in the domain containing the two common parts. 



5. Conclusion and discussion 

In this paper, we derived a number of probability distributions for the lengths and 
areas of triangles made of the three geodesies connecting three uniformly drawn random 
points, as well as of minimal separating loops. These laws are expected to be univer- 
sal features of the Brownian map and provide quantitative results characterizing the 
phenomenon of confiuence. This phenomenon is remarkable as it places the Brownian 
map half way between smooth surfaces and trees. In smooth surfaces, geodesies cannot 
merge and the three sides of a triangle only meet at their endpoints so that there are 
no common parts. On the contrary, in trees, the three sides of a triangle meet at a 
central common vertex so that there is no open part. As for a minimal separating loop 
on a tree, it corresponds generically to a back-and-forth travel to the above central 
common central vertex hence has no open part. On a smooth surface, depending on the 
shape of the surface, a minimal separating loop is either a back-and-forth travel along 
a geodesic, with no open part, or a simple curve with no common part. Having both 
open and common parts of non-zero length is a peculiarity of the Brownian map. 

Our approach consisted in obtaining discrete results for random quadrangulations 
and taking their scaling limit. So far we lack a general formalism which would allow to 



compute the same results directly in the continuum. Despite recent progress [ pT| , |22[] , the 
so-called Liouville field theory does not yet seem to be able to address such questions. 
Moreover, our results are restricted to the so-called universality class of pure gravity. 
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It would be desirable to extend them to other universality classes of random surfaces 
coupled to critical matter models (characterized by their central charge c, the pure 
gravity having c = 0) such as the celebrated Ising model (c = 1/2) |]2^. Discrete 



approaches based on bijections with blossom trees p^p6[| or labeled trees [p7| , p8[| exist 
for these problems but those could not be used so far to extract geometrical informations. 
Some of these models (with a central charge c > 1) are expected to behave like branched 
polymers, hence should have the geometry of trees described above. 



46 



References 



[1] V. Kazakov, Bilocal regularization of models of random surfaces, Phys. Lett. B150 
(1985) 282-284; F. David, Planar diagrams, two-dimensional lattice gravity and 
surface models, Nucl. Phys. B257 (1985) 45-58; J. Ambj0m, B. Durhuus and J. 
Frohlich, Diseases of triangulated random surface models and possible cures, Nucl. 
Phys. B257 (1985) 433-449; V. Kazakov, I. Rostov and A. Migdal Critical properties 
of randomly triangulated planar random surfaces, Phys. Lett. B157 (1985) 295-300. 

[2] for a review, see: P. Di Francesco, P. Ginsparg and J. Zinn-Justin, 2D Gravity and 
Random Matrices, Physics Reports 254 (1995) 1-131. 

[3] J. Ambj0rn, B. Durhuus and T. Jonsson, Quantum Geometry: A statistical field 
theory approach, Cambridge University Press, 1997. 

[4] J. F. Marckert and A. Mokkadem, Limit of normalized quadrangulations: the Brow- 
man map. Annals of Probability 34(6) (2006) 2144-2202, |arXiv:math.PR/0403398 



[5] J. F. Le Gall, The topological structure of scaling limits of large planar maps, invent. 

math. 169 (2007) 621-670, |arXiv:math.PR/0607567| . 
[6] J. F. Le Gall and F. Paulin, Scaling limits of bipartite planar maps are homeomor- 

phic to the 2-sphere, |arXiv:math.PR/0612315| . 
[7] G. Miermont, On the sphericity of scaling limits of random planar quadrangulations. 

Elect. Comm. Probab. 13 (2008) 248-257, |arXiv:0712.3687| [math.PR]. 
[8] J. Ambj0rn and Y. Watabiki, Scaling in quantum gravity, Nucl. Phys. B445 (1995) 

129-144. 

[9] J. Bouttier, P. Di Francesco and E. Guitter, Geodesic distance in planar graphs, 

Nucl. Phys. B663[FS] (2003) 535-567, |arXiv:cond-mat/ 03032 72| . 
[10] P. Chassaing and G. Schaeffer, Random Planar Lattices and Integrated SuperBrow- 

man Excursion, Probability Theory and Related Fields 128(2) (2004) 161-212, 

arXiv:math.CO/0205226| . 
[11] G. Miermont and M. Weill, Radius and profile of random planar maps with 

faces of arbitrary degrees. Electron. J. Probab. 13 (2008) 79-106, |arXiv:0706.3334 

[math.PR]. 

[12] J. Bouttier and E. Guitter, Statistics of geodesies in large quadrangulations, J. Phys. 

A: Math. Theor. 41 (2008) 145001 (30pp), |arXiv:0712.2T60| [math-ph]. 
[13] G. Miermont, Tessellations of random maps of arbitrary genus, |arXiv:0712.3688 

[math.PR] 



[14] J.-F. Le Gall, Geodesies in large planar maps and in the Brownian map, |arXiv:0804.3012 
[math.PR]. 

[15] J. Bouttier and E. Guitter, The three-point function of planar quadrangulations, J. 
Stat. Mech. (2008) P07020, |arXiv: 0805.2355] [math-ph]. 



47 



[16] M. Marcus and G. Schaeffer, Une bijection simple pour les cartes orientables, avail- 
able at pittp: / /www.lix.polytechnique.fr/Labo/Gilles.Scliaetfer/Biblio/l ; see also G. 
Schaeffer, Conjugaison d'arbres et cartes combinatoires aleatoires, PhD Thesis, Uni- 
versite Bordeaux I (1998) and G. Chapuy, M. Marcus and G. Schaeffer, A bijection 
for rooted maps on orientable surfaces, |arXiv:0712.3649| [math.CO]. 

[17] R. Cori and B. Vauquelin, Planar maps are well labeled trees, Canad. J. Math. 
33(5) (1981) 1023-1042. 

[18] J.-F. Delmas, Computation of moments for the length of the one dimensional ISE 
support, Elect. Journ. of Probab. 8(17) (2003) 1-15. 

[19] M. Bousquet-Melou, Limit laws for embedded trees. Applications to the integrated 
superBrownian excursion, Random Structures and Algorithms 29(4) (2006) 475- 
523, |arXiv:math/0501266| [math.CO]. 

[20] J. Bouttier, Physique statistique des surfaces aleatoires et combinatoire bijective 
des cartes planaires, PhD Thesis (2005). 



[21] B. Duplantier and S. Sheffield, Liouville Quantum Gravity and irPZ, |arXiv:0808.1560 
[math.PR]. 

[22] F. David and M. Bauer, Another derivation of the geometrical KPZ relations, 

^rXiv:08l0.^83B| [math-ph]. 
[23] D. Boulatov and V. Kazakov, The Ising model on a random planar lattice: the 

structure of the phase transition and the exact critical exponents, Phys. Lett. B186 

(1987) 379-384. 

[24] M. Bousquet-Melou and G. Schaeffer, The degree distribution in bipartite planar 

maps: application to the Ising model, |arXiv:math.CQ/0211070| . 
[25] J. Bouttier, P. Di Francesco and E. Guitter. Combinatorics of Hard Particles on 

Planar Graphs, Nucl.Phys. B655 (2003) 313-341, |arXiv:cond-mat/021116'g . 
[26] J. Bouttier, P. Di Francesco and E. Guitter. Combinatorics of bicubic maps with 

hard particles, J.Phys. A: Math.Gen. 38 (2005) 4529-4560, |arXiv:math.CO/0501344 . 
[27] J. Bouttier, P. Di Francesco and E. Guitter. Planar maps as labeled mobiles, Elec. 

Jour, of Combinatorics 11 (2004) R69, |arXiv:math.CO/0405099 . 
[28] J. Bouttier, P. Di Francesco and E. Guitter. Blocked edges on Eulerian maps and 

mobiles: Application to spanning trees, hard particles and the Ising model, J. Phys. 

A: Math. Theor. 40 (2007) 7411-7440, |arXiv:math.CU/0702097i 



48 



